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*5 '. Abstract 

The semi-exponential basis set of radial functions (A.M. Frolov, Physics Letters A 374, 
2361 (2010)) is used for variational computations of bound states in three-electron atomic 

■ systems. It appears that semi-exponential basis set has a substantially greater poten- 
tial for accurate variational computations of bound states in three-electron atomic systems 

>> 

' than it was originally anticipated. In particular, the 40-term Larson's wave function im- 

proved with the use of semi-exponential radial basis functions now produces the total energy 
-7.4780581457 a.u. for the ground l 2 ^— state in the °°Li atom (only one spin function \i = 
af3a — f3aa was used in these calculations). This variational energy is very close to the exact 

in 

CN ' ground state energy of the °°Li atom and it substantially lower than the total energy obtained 

OO 

■ with the original Larson's 40-term wave function (-7.477944869 a.u.). 
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I. INTRODUCTION 



In this study we perform variational calculations of bound states in three-electron atomic 
systems. A basis set of semi-exponential radial functions [1] is extensively used in our 
calculations. The main goal is to solve the non-relativistic Schrodinger equation H^> = E^fr, 
where E < and bound state wave function \1/ has the unit norm. The general non- 
relativistic Hamiltonian H of the three-electron atomic problem is (see, e.g., [2|) 



H 



Ti 2 r _o _o _o m P _oi Qe 2 Qe 2 Qe 2 e 2 e 2 e 2 



2m P 



vi + vi + v 2 3 + -^vj - - — - — —+—+—+— (i) 

M J ri4 r 24 r 34 r V2 r i3 r 23 

where Ti — is the reduced Planck constant, m e is the electron mass and e is the electric 
charge of electron. In this equation and everywhere below in this study the subscripts 
1, 2, 3 designate the three electrons e~~, while the subscript 4 denotes the heavy nucleus 
with the mass M (M 3> m e ) and positive electric (nuclear) charge Qe. The notations 
rij =| r i — Tj |= Tji stand for the six interparticle distances (= relative coordinates) defined 
in an arbitrary four-body system and (i — 1, 2, 3, 4) are the Cartesian coordinates of the 
four point particles. In Eq.([T]) and everywhere below in this work we shall assume that (ij) 
= (ji) = (12), (13), (14), (23), (24), (34). Below only atomic units h = 1, | e |= 1, m e = 1 
are employed. In these units the explicit form of the Hamiltonian H, Eq.((T]), is significantly 
simplified. 

The main attention in this work is focused on numerical calculations of the ground (dou- 
blet) 1 2 S(L = 0)— state (or l 2 Si(L = 0)— state) of the three-electron Li atom with the 
infinitely heavy nucleus, i.e. the °°Li atom. Considerations of other three-electron atoms, 
ions and various positron containing atomic systems (e.g., HPs) can be performed absolutely 
analogously (for more detail, see [lj) and, therefore, these systems will not be considered 
here. 

The problem of highly accurate calculations of the bound states in three-electron atomic 
systems has attracted continuing attention. The first calculations of the Li atom with 
the truly correlated wave functions were performed in 1936 [3j. A brief reviews of such 
calculations can be found in jj] (earlier works) and j^J (references up to 1997 are mentioned). 
The current bibliography on this subject includes almost one thousand references and is 
increasing rapidly. The classical Hylleraas method (see, e.g., [5]) is the method capable of 
predicting the most highly accurate wave functions for bound states of three-electron atomic 
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systems. In this method to produce the most highly accurate wave functions (e.g., for the 
ground (doublet) l 2 S'-state in the Li atom) one needs to use many thousands of Hylleraas 
basis functions. The use of extremely large basis sets is very inconvenient in many actual 
cases, since it produces a number of computational problems. It is clear that the classical 
Hylleraas method cannot be used to construct both compact and accurate variational wave 
functions for three-electron systems. Indeed, it contains essentially no control parameters 
which can be optimized by increasing the overall efficiency of the method. 

An alternative approach to variational bound state calculations in three-electron atomic 

n 

systems was proposed in This approach is based on the use of semi-exponential varia- 
tional wave functions l| and it allows one to construct very compact and accurate variational 
wave functions for arbitrary three-electron atomic system. Each of the semi-exponential 
basis functions depend upon all six interparticle coordinates ri 2 , ri 3 , . . . , r 34 A high 
efficiency of this new approach in actual applications and its superiority over the classi- 
cal Hylleraas expansion was demonstrated in jl]. As follows from the results of this study 
the semi-exponential basis set has a substantially greater potential for highly accurate varia- 
tional computations of bound states in three-electron atomic systems than it was anticipated 
earlier 

II. VARIATIONAL WAVE FUNCTION 

The variational wave function of the doublet S(L = 0)— states of the three-electron Li 
atom is written in the following general form 

^l=o = ifiL=o{A; {rij})(a/3a - f3aa) + (p L=0 (B; {r ij }){2aa(5 - /3aa - a/3a) (2) 

where -0£ =O (A; {r^}) and (pL=o(B; {r^}) are the two independent spatial parts (= radial 
parts) of the total wave function. Each of these two radial functions is, in fact, a radial 
factor (for states with L = 0) in front of the corresponding three-electron spin functions 
Xi = a/3a — (3aa and \2 = 2aa/3 — (3aa — a/3a. Here the notations a and /3 are the one- 
electron spin- up and spin-down functions, respectively (see, e.g., 6]). The notations A and B 
in Eq. (T5]) mean that the two sets of non-linear parameters associated with radial functions 
and cj) can be optimized independently. In the general case, each of the radial basis functions 
explicitly depends upon all six interparticle (relative) coordinates r 12 , r 13 , r 2 3, r 14 , r 2 4, r 34 . It 
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is clear that in actual variational calculations only one spin function, e.g., the \i function, 
need be used. Another useful trick (so-called 'doubling') is based on the use of the same set 
of non-linear parameters in the two radial parts in Eq.(J2]). 



In our earlier work [l| we have introduced an advanced set of radial basis functions for 
three-electron atomic calculations. In |l| this set was called the semi-exponential basis set. 
In general, the semi-exponential variational expansion of the radial function if)L = o(A; {ry}) 
is written in the form 

N 



ij L=0 (A; { rij }) = £ C fe r 2 "3 l(fc) r- (fc) rr 2 3(fc) rr 4 l(/£) r 2 7 (fc) r 3 7 (fe) e^(-a fc r 14 - /3 fc r 24 - 7fe r 34 ) (3) 

k=l 



where akihilk (k = 1,2, . . . , N) are the varied non-linear parameters. The presence of 
the varied non-linear parameters in Eq.flBJ) is the main and very important difference with 
the traditional Hylleraas variational expansion (see, e.g., (7|) for which in Eq.flH]) we always 
have ax = ... = &n, Pi = ■ ■ ■ = 0n and 71 = . . . = 7^. Note that all matrix elements of 
the Hamiltonian, Eq.flT]), and overlap matrix needed in computations with the use of the 
semi-exponential basis, Eq.(j3]), contain the same three-electron integrals which arise for the 
usual Hylleraas expansion (for more detail, see In other words, numerical calculation 
of all matrix elements with semi-exponential functions is no more difficult problem, than 
for the traditional Hylleraas radial functions. This also simplifies numerical computation 
of the bound state properties (i.e. expectation values) in the semi-exponential basis set. 
Our algorithms used in calculations of all required matrix elements is based on the Perkins 
formula for three-electron integrals [8] in relative coordinates. Note also that all calculations 
in this work have been performed with the use of standard quadruple precision accuracy (30 
decimal digits per computer word). 

In actual atomic systems any many-electron wave function must be completely an- 
tisymmetric upon all electron variables, i.e. upon all electron spatial and spin vari- 
ables. For three-electron atomic wave function this requirement is written in the form 
^4i23*(l,2,3) = -#(1,2,3), where $ is given by Eq.© and A e is the three-particle (= 
electron) antisymmetrizer A e = e — A2 — A3 — A3 + A23 + A32- Here e is the identity 
permutation, while A? is the permutation of the z-th and j-th particles. Analogously, the 
operator Ajfc is the permutation of the i-th, j-th and k-th particles. In actual computations 
antisymmetrization of the total wave function is reduced to the proper antisymmetrization 



of corresponding matrix elements (for more detail, see, e.g., jlj). Each of these matrix ele- 



merits is written in the form | O \ ^f), where O is an arbitrary spin-independent quantum 
operator which is truly symmetric upon all interparticle permutations. The wave function 
Eq.([2]), contains the two different radial parts ip and <f>. By performing the integration 
over all spin coordinates from here one finds the four spatial projectors V^,V^<j, = and 
V^rj, presented in In fact, the explicit form of the and projectors given in [l| 
must be corrected (there is an obvious misprint in the formulas given in 

= \ (A 3 - A 3 + A 23 - A 32 ) (4) 

'P(/>i> — 2 (-^13 — A3 + A23 — A32) (5) 

For an arbitrary truly symmetric spin-independent operator O each of these four projectors 
produces matrix elements \ O \ ty) of the correct permutation symmetry (for doublet 
states) between all three electrons. The explicit formulas for all matrix elements obtained 
with the radial basis functions, Eq.(j3]), and for three-electron integrals needed in calculations 
can be found in 



III. CALCULATIONS 



Let us apply the semi-exponential variational expansion, Eq.(|3]), to numerical calcula- 
tions of the ground l^-state in the three-electron °°Li atom. In this study we consider the 
two variational wave functions: (a) the wave function which contains 28 radial basis func- 
tions, Eq.(j3]), and (b) the wave function which includes 40 radial basis functions, Eq.(j3]). 
The results (in atomic units) obtained with these two trial wave functions can be found in 
Table I. Tables II and III contain the corresponding radial basis functions, Eq.(|3]), i.e. the 
powers ni(k), ^(fc), n^k), mi(k), rri2{k), m^{k) of six radial variables r^, ri3, r23, r^, r2A, 
and optimized non-linear parameters ctk,Pk,lk- As follows from Table I our variational en- 
ergies obtained for the ground l 2 S'-state in the °°Li atom with the use of semi-exponential 
variational expansion, Eq.([2]), are substantially lower than the corresponding energies deter- 
mined for this state with the same Hylleraas wave function [7[ . Note also that the non-linear 
parameters used in our method (in Eq.(j3])) are constantly varied. Therefore, it is hard to 
say that the total energies obtained in some calculations are 'final'. Formally, based on the 
known convergence rate(s) for our data and by using a few extrapolation procedures we 
can approximately evaluate the limits to which our variational energies will converge, if we 
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could perform an infinite number of variations for the non-linear parameters in Eq.(j3]). Such 
limits for the total energies are shown in the fourth column of Table I. These values indicate 
that, e.g., our 40-term variational wave function can produce, in principle, very accurate 
variational energies, if the optimization of non-linear parameters in Eq.flH]) will continue. 

In {2I Larsson proposed a simple (but useful!) trick which allows one to increase the 
overall accuracy of the trial (doublet) wave function. Later this trick was called 'doubling' 
of the wave function and it was used practically in all calculations of the bound doublet 
states in three-electron atomic systems. The idea of doubling is simple and transparent. 
If we already know the radial function constructed for one spin configuration, e.g., for 
Xi = oi(5a — (3aa from Eq.(l2"l), then we can use exactly the same radial basis function for 
another spin configuration \2 — 2ctaf3 — fiaa — aj3a. Formally, it doubles the total number 
of basis functions in the trial wave function. According to the variational principle the 
total variational energy can only decrease during such a procedure. The problem of linear 
dependence of basis functions is avoided in this procedure, since the two spin functions 
Xi = a(5a — (5aa and X2 = 2aa/3 — (3aa — a(3a are independent of each other. In fact, 
for Hylleraas variational expansion the 'doubling' does not work properly, since there are 
obvious linear dependencies between different radial basis functions in those cases when some 
non-linear parameters coincide with each other (for more detail, see [?j]). In semi-exponential 
variational expansion Eq.(j3]) all optimized non-linear parameters are independent of each 
other. Therefore, the coincidence of the pre-exponential factors in Eq.([3]) is not crucial and 
does not mean that such basis functions are linearly dependent. This drastically simplifies 
the actual 'doubling' for Eq.(|3]). The energies obtained with the use of 'doubling' of our 
variational wave functions can be found in the fifth column of Table I. It is clear that the 
non- linear parameters (or parameters in the exponents in Eq.Q) from the second part of 
the total wave function are not optimal, i.e. for all terms which contain basis functions 
with numbers i > 41 the non-linear parameters are not optimal. These 120 (= 40 x3) non- 
linear parameters can be re-optimized and this drastically improves the overall quality of 
the total wave function. For instance, approximate re-optimization of the last 40 non-linear 
parameters in the wave function gives the ground state energy -7.4780583419 a.u., which is 
much better than the 'doubling' energy (-7.4780581691 a.u.) from Table I. Note also that 
our value of the total energy of the ground state in the Li atom is better than the total 
energy obtained in [10] with the use of 352 basis functions. 
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As follows from Table I the doubling is not an effective approach for our trial wave 
functions with the carefully optimized non-linear parameters. However, we can modify the 
original idea of doubling into something new which is substantially more effective in actual 
computations. To illustrate one of such modifications, let us assume that we have con- 
structed a 40-term variational wave function, Eq.flH]) which contains 40 x 3 = 120 carefully 
optimized non-linear parameters a%, fix, 71, . . . , a^o, P40, 740- At the second step of our pro- 
cedure we can add forty additional basis functions with the same pre-exponential factors 
^23 ^13 r^l^r^^V^^V^ 3 ^, but slightly different exponents in Eq.Q. In reality, these 
new exponents have been chosen quasi-randomly from three different intervals, e.g., 

a i+m = a i + 0.0057 + 



A+40 = A + 0.0063 



2 

i{i+l)V3 
2 



7i+ 4o = 7 ? + 0.0049 • + )) 



where i = 1, 2, . . ., 40 and \\ x )) designates the fractional part of the real number x. Small 
deviations of these new exponents from the known 'optimal' values (i.e. from 0^,^,7$, 
where 1 < i < 40) produce the extended wave function of 'almost optimal' quality. On 
the other hand, even these small differences between exponents allows one to avoid a linear 
dependence between basis vectors in Eq.Q. Obviously, this procedure can be repeated a 
number of times. This allows one to construct very accurate trial wave functions which 
contain not only 80, but 400, 800 and even 2000 basis functions with almost 'optimal' non- 
linear parameters. 

It is very interesting to perform variational calculations of the ground state of the °°Li 



atom with the use of 60-term wave function constructed from the analo 



erous 60-term wave 



7] with that wave 



Larsson's wave function The variational total energy obtained in 
function and one spin function xi was -7.4780103597 a.u. Our 60-term trial wave function 
with one spin function Xi constructed in [I] from the same radial basis functions corresponds 
to the substantially lower total energy E = -7.478057561 a.u. The current total energy 
is -7.4780597045 a.u. (only one spin function xi is used in our calculations). Note that 
our current total energy rapidly decreases with almost constant rate ~ 0.5 • 10~ 7 a.u. per 
optimization cycle, i.e. per one variation of all 180 (= 60 x 3) non-linear parameters in 



the trial wave function. This total energy is better than the values obtained in 



llj with 



the use of 503 Hylleraas basis functions (selected radial configurations were used in such 
calculations). Moreover, this total energy is slightly better than the value obtained in 12]. 

The doubling of our 60-term wave function produces the total energy -7.4780597761 a.u. 
This energy is close to the known 'exact' answer -7.478060323904 a.u. [9|. We expect that 
after an infinite number of variations of the non-linear parameters the total energy of the 
°°Li atom obtained with our 60-term radial function and one spin function will converge to 
the value -7.4780603(3) a.u. which is very close to the actual ground state energy. It will be 
an outstanding result to obtain the value lower than -7.4780602 a.u. for the total energy of 
the ground l 2 ^— state in the °°Li atom by using only 60-term variational wave function. 



IV. DISCUSSIONS AND CONCLUSION 

The semi-exponential variational expansion [l| is applied for bound state calculations 
of three-electron atomic systems. It is shown that this variational expansion allows one to 
construct compact and accurate variational wave functions for three-electron atomic systems. 
Currently, the use of semi-exponential radial basis functions is the best way to produce 
compact and accurate wave functions for three-electron atomic systems. The total energies 
obtained in this study for the ground l 2 S'-state of the °°Li atom are more accurate than 
our earlier results from 11 and substantially more accurate than the original Larsson's wave 

n 

function [7J with the same number of terms. 

The results indicate clearly that our semi-exponential variational expansion, Eq.flHJ, has 
a substantially greater potential for variational bound state calculations in three-electron 

n 

atomic systems than we have anticipated originally Currently, we continue the process of 
numerical optimization of the non-linear parameters in our trial wave functions constructed 
with the use of semi-exponential variational expansions. Note that the choice of Larsson's 
wave function(s) as the first approximation to the semi-exponential variational expansion is 
not crucial for our method. Many other choices are also possible. For instance, our next 
step will be re-optimization of the non-linear parameters in the 352-term wave function used 



in 



10] . It will take some time, but with such a wave function we hope to produce the total 



energy which is close to the known 'exact' energy for the ground state of the °°Li atom (E 
= -7.478060323904 a.u.). 

When our research of the semi-exponential variational expansion started we could not 



expect such impressive results. Currently, the variational wave functions constructed with 
the use of semi-exponential variational expansion are of great interest for many scientific 
problems which include the Li atom and other atomic and quasi-atomic three-electron sys- 
tems. In particular, we have made numerous improvements in our original computer code 
Our next goal is to generalize the semi-exponential basis to the four-electron atomic 
problems. The semi-exponential variational expansion for four-electron atomic systems can 
be used to obtain compact and accurate wave functions of the ground singlet state and 
triplet 2 3 5'— states of the Be atom and Be-like ions. It is clear that our method can also be 
used for rotationally excited states in atomic systems, i.e. for states with L > 1, where L is 
the electron angular momentum. 

In conclusion, it should be mentioned that numerous attempts to improve the overall 
quality of radial basis functions for three-electron atomic systems (in general, for four- 
body Coulomb systems) started almost 20 years ago (see, e.g., [ijj], [U^, [3]). In part, it 
was a reaction to low efficiency of the traditional Hylleraas variational expansion for such 
systems. Indeed, the total number of terms in modern versions of Hylleraas variational 
expansions already exceeded 13,000 and still growing. For the Li atom (ground state) such 
wave functions allows one to determine 12 - 13 stable decimal digits in the total energy. 



Comparison of a large number o: 
can be found in Table II from 



numerical results for the ground state of the lithium atom 



131 ] . The new Hylleraas wave function with 65,000 terms 



will allow one to obtain « 15 stable decimal digits in the total energy. It is clear that 
such a method cannot be considered as a reasonable and appropriate approach for accurate 
variational computations of the bound states in four-body systems. 

An alternative approach for construction of highly accurate four-body wave functions was 



proposed in 



14j . 15j, 16|. It is based on the closed analytical formulas derived in [l7] for 



the basic four-body integral. As it follows from actual calculations this method allows one 
to obtain very accurate bound state energies essentially for all Coulomb four-body systems, 
including °°Li atom (E = -7.47806025114 a.u., 50 radial basis functions (exponents) with the 
spin function xi phis 19 radial basis functions with the spin function X2 (their definitions are 
given above), bi-positronium Ps2 (E = -0.51600377267 a.u., 50 exponential basis functions) 
and °°HPs (E = -0.78914861151 a.u., 50 exponential basis functions). Very good values 
for the total (bound) state energies have also been obtained for the ppufi and dd\i\i four- 
body bi-muon systems. These total energies can be improved even further by using more 
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careful optimization of the non-linear parameters in the trial wave functions. For accurate 
computations of the H2 molecule and other similar (molecular) systems the method [ljj] 
must be modified to include the complex values for some non-linear parameters (as was 

n 

done in |16|). 

Note, however, that relations between these variational expansions based on the formula 
from [l7| and traditional Hylleraas variational expansion are extremely complicated. Further 
analysis shows that the source of these difficulties is directly related with the fundamental 
properties of the four-body perimetric coordinates fis|. rather than with the lack of 'good 
representations' for the formula Eqs.(2.1) - (2.9) from [l7|. As follows from the properties of 
the four-body perimetric coordinates [18j it is very hard to represent any of these compact 
and accurate 'exponential' wave functions in terms of the Hylleraas basis functions. The 
overall complexity of problems arising here is comparable with the difficulties which can 
be found in the original problem. Many advantages of the 'exponential' wave functions, 
e.g., their very accuracy, can be lost during this procedure. In many cases, it is simpler 
to re-calculate the corresponding energies from the very beginning by using the traditional 
Hylleraas variational expansion. This is the main reason why we do not want to discuss here 
some recent developments in the four-body exponential and related variational expansions 
(see, e.g., [19|). 



Our approach allows one to produce compact and accurate wave functions which are 
easily related with the Hylleraas basis set of radial functions. This means that all such 
compact and accurate wave functions can directly be used in computations of various bound 
state properties and transition probabilities. The second advantage of the semi-exponential 
variational expansion follows from relatively simple formulas for all matrix elements needed 
in highly accurate computations and optimization processes. 
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TABLE I: The total energies E (in atomic units) of the 1 2 S(L = 0)— state in the °°Li atom. N 
designates the number of basis functions used. 

N £(Ref.[4]) £(Eq.(3)) £ a (Eq.(3)) £(Eq.(3); doubling) 
28 -7.477885105 -7.4780363801 -7.4780368(3) -7.4780365786 
40 -7.477944869 -7.4780579457 -7.4780595(5) -7.4780580161 
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TABLE II: An example of the trial wave function constructed with the use of N = 28 semi- 
exponential radial basis functions. This wave function produces the total energy E = -7.4780363801 
a.u. for the ground l 2 ^— state of the °°Li atom. Only one electron spin-function \i = a P a ~ P aa 
was used in these calculations. 
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0.672095727321017E+00 


16 








4 








1 


0.113723992478698E+02 


0.508008324231496E+01 


0.636162427376743E+00 


17 











2 


2 


1 


0.368595304619283E+01 


0.278359204737422E+01 


0.620935860917418E+00 


18 








1 


1 


1 


1 


0.357466607750891E+01 


0.343630387369804E+01 


0.698563354607838E+00 


19 








1 


2 





1 


0.365095130448352E+01 


0.317945977760753E+01 


0.652269908234656E+00 


20 








1 


3 





1 


0.290005170216145E+01 


0.289825235917765E+01 


0.685753232846547E+00 


21 


2 

















0.280340788922781E+01 


0.225117021991187E+01 


0.990926897297406E+00 


22 


1 


1 














0.275272575017714E+01 


0.276113990248958E+01 


0.802457922378512E+00 


23 


1 











2 





0.606571100003296E+01 


0.326870035708320E+01 


0.453868977046091E+00 


24 


1 








1 


1 





0.293743657648302E+01 


0.226578275529031E+01 


0.669408905271258E+00 


25 

















4 


0.502952034298549E+01 


0.276772288773810E+01 


0.950599883558449E+00 


26 








1 











0.868992405342042E+01 


0.374899612568756E+01 


0.516228416453001E+00 


27 








1 








2 


0.433237156102395E+01 


0.278286578531691E+01 


0.112727477282094E+01 


28 











1 








0.221163211952812E+01 


0.498938978770980E+01 


0.647628325318278E+00 



13 



TABLE III: An example of the trial wave function constructed with the use of N = 40 
semi-exponential radial basis functions. This wave function produces the total energy E = - 
7.47805542591 a.u. for the ground l 2 ^— state of the °°Li atom. Only one electron spin-function 
Xi = af3a — j3aa was used in these calculations. 



N 


"i 


»i2 n 3 mi m-2 rri3 


a 


P 


7 


1 

















1 


0.310888325274582E+01 


0.290883960347876E+01 


0.845418768108028E+00 


2 











1 





1 


0.196342255489068E+01 


0.347479833068431E+01 


0.163523383561211E+01 


3 











1 


1 


1 


0.278598083293758E+01 


0.282092658543858E+01 


0.679130247211775E+00 


4 











2 





1 


0.287878476756554E+01 


0.284028201233856E+01 


0.652121055380736E+00 


5 








1 








1 


0.296290423881450E+01 


0.281951549636015E+01 


0.660626319649503E+00 


6 








2 








1 


0.356029841316233E+01 


0.298528920333113E+01 


0.672210253693381E+00 


7 




















0.281820369226934E+01 


0.307242548167154E+01 


0.614180061849861E+00 


8 


1 

















0.129540716567888E+01 


0.388006326151676E+01 


0.588477295803696E+00 


9 

















2 


0.312600561158126E+01 


0.298186247192942E+01 


0.711199626874967E+00 


10 


1 











1 





0.309479768715561E+01 


0.205912634163624E+01 


0.692766055505988E+00 


11 








3 








1 


0.421833159497161E+01 


0.323537515175651E+01 


0.683447854994486E+00 


12 


1 














1 


0.317935937791374E+01 


0.246672317399557E+01 


0.128821985777175E+01 


13 

















3 


0.206773619817191E+01 


0.258067622963825E+01 


0.102117374065975E+01 


14 








1 


1 





1 


0.407929602782570E+01 


0.308726834443555E+01 


0.145857028906218E+01 


15 











3 





1 


0.290479112486719E+01 


0.342131088097724E+01 


0.654819087759473E+00 


16 








4 








1 


0.759200717376427E+01 


0.403600107010476E+01 


0.654094476238819E+00 


17 











2 


2 


1 


0.342043341882141E+01 


0.300593652535184E+01 


. 75898039071 1669E+00 


18 








1 


1 


1 


1 


0.337729059232100E+01 


0.336098394507054E+01 


0.689553038540639E+00 


19 








1 


2 





1 


0.376783902187767E+01 


0.308351505579080E+01 


0.693591996711379E+00 


20 








1 


3 





1 


0.341103485088192E+01 


0.390874343351408E+01 


0.353215493102457E+01 


21 


2 

















0.187375370334190E+01 


0.292204564107203E+01 


0.134200544832659E+01 


22 


1 


1 














0.203602651426972E+01 


0.251346948989043E+01 


0.139512883704298E+01 


23 


1 











2 





0.307587342243077E+01 


0.289804349780859E+01 


0.799088792948691E+00 


24 


1 








1 


1 





0.237469500100334E+01 


0.411142232808047E+01 


0.113515551805206E+01 


25 

















4 


0.306630970495110E+01 


0.358372314305121E+01 


0.865024537355183E+00 


26 








1 











0.652799588882586E+01 


0.391736310812260E+01 


0.504690239191085E+00 


27 








1 








2 


0.523338852365502E+01 


0.400156315581741E+01 


0.119987716494855E+01 


28 











1 








0.303300703612990E+01 


0.269004008151858E+01 


0.132439178382705E+01 


29 











1 





2 


0.347252909237306E+01 


0.351142344973601E+01 


0.841746594395407E+00 


30 


1 





1 











0.246738217080682E+01 


0.268143054665150E+01 


0.930487658749999E+00 


31 


2 











1 





0.313058760564443E+01 


0.295867936877980E+01 


0.903526045564038E+00 


32 


1 











1 


1 


0.324981907924628E+01 


0.322474651197190E+01 


0.896978954638563E+00 


33 


3 

















0.333729310041323E+01 


0.277869134235698E+01 


0.100059195659526E+01 


34 


2 














1 


0.345708067187951E+01 


0.277492748358036E+01 


0.104744023031760E+01 


35 








5 








1 


0.621080396345153E+01 


0.394717905188376E+01 


0.771357719659273E+00 


36 











4 





1 


0.617132034807857E+01 


0.257389911129776E+01 


0.405076099224201E+01 


37 








1 


4 





1 


0.384658211093920E+01 


0.426365221909316E+01 


0.963468025067806E+00 


38 











5 





1 


0.358540307325767E+01 


0.426395800491431E+01 


0.120557474663118E+01 


39 








2 


1 





1 


0.290948705798909E+01 


0.342780417090815E+01 


0.803444693562464E+00 


40 








2 


2 





1 


0.348532012166237E+01 


0.473382479410345E+01 


0.871338172104194E+00 



14 



